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ABSTRACT 


The  solution  is  presented  for  the  problems  of  edge 

dislocation  in  hexagonal  crystals  with  long  range  , 

interatomic  interactions.  The  field  equations  of  i 

nonlocal  elastic  solids  are  employed  to  determine  i 

the  stress  fields  and  the  elastic  energy  for  an 
edge  dislocation  in  the  basal  plane.  Classical 
stress  and  energy  singvilarities  are  found  not  present 

in  the  nonlocal  solutions.  Stress  distribution  , 

is  calculated  and  maximum  shear  stress  is  given  for  I 

various  hexagonal  materials.  Theoretical  shear  i 

stress  to  initiate  a dislocation  having  a Burger's  * 

vector  of  one  atomic  distance  is  calculated  and  i 

found  to  be  in  the  acceptable  range  known  from  the  i 

lattice  dynamic  calculations. 


1.  INTRODUCTION 

It  is  well-known  that  the  classical  elasticity  solution  of  the 
edge  dislocation  contains  stress  and  energy  singularities  in  the 
"core  region",  cf.  [l].  In  several  previous  papers  (e.g.,  [2], 

[3],  [^])  we  have  shown  that  the  solutions  of  various  Volterra 
dislocations  based  on  the  nonlocal  elasticity  theory  do  not  contain 


« 
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those  singularities.  This  recent,  theory  [^,6]  models  the  elastic 
materials  much  more  satisfactorily  in  tiiat  the  effect  of  long  range 
interatomic  interactions  are  taken  into  account.  It.  seems  that  no 
artifice  such  as  introducing  various  atomistic  models  to  estimate 
the  stress  and  energy  in  t)ie  core  region  is  necessary.  Moreover,  as 
a continuum  theory  all  problems  can  be  reduced  to  I'oundary-initial 
value  oroblems. 

The  discussion  of  the  dislocation  problems  in  anisotropic 
solids  is  not  a trivial  extension  even  in  its  classical  frame  of 
roCprence.  Moreover  because  of  the  orientational  effects  the  state 
of  stress  and  elastic  energy  arc  affected  considerably.  Constapient ly 
the  criteria  for  failure  or  the  generation  of  dislocations  need  new 
investigations.  The  raison  d'etre  of  the  present  naper  stems  from 
these  considerations.  In  Sect.ion  ? we  present  a brief  summary  of 
the  field  equations  of  the  nonlocal  elasticity  tlieory.  Tn  Section 
3 we  obtain  the  solution  of  the  edge  dislocation  problem  leading 
to  stress  and  energy  fields  and  in  Section  we  specialise  these  results 
to  the  isotropic  crystals.  In  Section  5,  some  results  of  comn\jter 
calculations  are  presented  and  maximum  stresses  that  cause  a single  edge 
dislocation  in  several  hexagonal  crystals  (Mg.  Apatite,  Cd,  ”n ) are 
calculated,  '.fhe  distribution  of  normal  and  shear  stresses  along  radial 
line  r,  6=0,  and  as  a function  of  the  polar  angle  0 for  a fixed  r are 
calculated,  flince  no  stress  and  energy  sing\ilarity  occ\ir,  the  naxin\im 
stress  hypothesis  may  be  used  to  calculate  the  theoretical  (cohesive^ 
stress.  The  results  are  in  the  accepted  range  known  from  tlie  atomic 


theory. 
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2.  FORMULATION 


In  previous  papers  [2,i*]  we  have  shown  that,  under  some  very 
general  conditions,  the  solution  of  elastostatic  probleias  in  linear 
nonlocal  elasticity  can  be  reduced  to  the  solution  of  the  classical 
Navier's  equation,  however  the  stress  field  is  calculated  by 


(2.1) 


j a(x'-x)aj^^(x' )dv(x' ) 
V 


where  is  given  by  the  classical  Hooke's  law  which  for  the 
hexagonal  crystals  can  be  arranged  into  the  form,  cf.  Fig.  la. 
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Here  is  the  linear  strain  tensor  given  by 


where  u^  Since  for  the  hexagonal  crystals  , 

the  number  of  independent  elastic  constants  are  five. 


k 


The  attenuation  function  a{x'-x)  sug^^ested  in  our  previous  work 
is  of  the  form 

(2. It)  a = Cl  exp[-(k  /a)^(x’-x.  )(x'-x.  )-(ko/a)^(x^-x^)^] 

where  B is  summed  over  B=1  and  B=3.  Here  a Is  the  lattice  parameter 

k^  and  k^  are  two  consteints  which  govern  the  range  of  attenuation  of 

the  Interatomic  attractions.  The  constant  a is  determined  from  the 

o 

normalization  condition  on  a; 

(2.5) 

so  that 

(2.6)  a = kfk.,/it^^^  b}  . 

Now  a problem  in  nonlocal  elasticity  is  reduced  to  determining  the 
displacement  field  '^(x)  by  solving  the  Navler's  equation  obtained 
by  combining  (2.2)  and  (2.3)  with  the  Cauchy's  equation 

’M.k'"  ' '<,..1.2,3 

where  as  usual  repeated  indices  are  summed  over  (1,2,3).  Once  ^(x) 
is  determined  one  can  then  calculate  strain  from  (2.3)  and  the  stress 
field  tj^i^  by  using  (2.2)  in  (2.1).  We  now  apply  this  program  to  solve 
the  problem  of  the  edge  dislocations. 


J a(x'-x)  dv(x' ) = 1 
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edgf:  dislocation 

The  straight  edge  dislocation  in  a hexagonal  crystal  la  possible 

in  the  basal  plane  y=const.  in  the  x-directlon  Fig.  la.  We  vision 

such  a dislocation  by  cutting  a cylinder  along  a radial  plane  and 

pulling  lower  surface  by  a constant  amount  b (called  Burger's  vector) 

relative  to  the  upper  surface  and  welding  the  two  surfaces.  Fig.  lb. 

The  solution  of  this  problem  in  classical  elasticity  is  well-known 

(cf.  [l],  p.  U2C ) . The  displacement  field  (u  ,u  ,0)  satisfving  Navier's 

X y ■ 

equation,  is  given  by 
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all  other  components  of  and  u^  vanish.  Hero  we  used  (x,y,?.) 
for  subscripts  (1,2,3)  and  set 


6 


I li 


(3.3) 
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u H X/(C  sin24')  , = (C^^-C^^)cos4. 


Introducing  the  cylindrical  coordinates  by  (Fig- 


(3.U)  X = r cose  , y = r sinO  , z = z 


r’  cose'  , y'  = r'  sine'  , z*  = z 
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U.sitip  tlie  relations 
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(3.8)  nay  be  written  as 
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The  attenuation  funotion  n(x'-x)  in  cylindrical  coordinates  ncquirop 
the  form 


(3.11)  a(x*-x)  = IT  ^^‘■(k!^k^/a^)exi‘'[-(k^/n‘‘^)(s'-n)''l 
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We  now  substitute  (3.1*)  into  (3.2),  the  result  and  (3.1l)  into  (3.10). 
This,  for  each  of  the  stress  components,  leads  to  a triple  integral 
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over  the  domain  (0^r'<“’,  0^(''<2tt,  -®<2'<oo).  The  integrations  over  r' 
and  z'  are  tedious  but  can  be  carried  out,  however  the  integration  over 
0'  will  have  to  be  done  numerically.  Leaving  the  details  of  these 
calculations  we  give  the  results: 
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The  total  strain  energ\'  is  piven  by 
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where  L is  the  length  of  the  cylinder  and  R is  the  radius.  Substituting 

for  e , e.,^  and  e „ calculated  from  (2.3),  in  the  same  way  as  in  the 
rr  99  r9 

case  for  , we  write 
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where 
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h . ISOTROPIC  SOLID 

In  the  case  of  isotropic  solids  we  have 


(U.l)  M = W/(l-v)  , K=i  , ^=^<12=2  , q3=‘*  , qi,=o  , 
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where  C ts  the  Killer's  constant  and 


(It.l*)  p = kr/a  , — ‘11  •«  = 0.''T7:'’. . . 


Tliese  results  are  In  complete  apreement  with  t.hose  plven  In  (ll. 


s.  nrscainTON 

Components  of  the  stress  tensor  are  pjott.ed  in  Kip.  1 as  a funct  ion 

of  the  polar  angle  0 for  a fixed  radial  distance.  Tt  ts  observed  that 

t and  t„-  have  the  same  shape  with  their  extrema  oceurrlnp  at  tdie  siune 
rr  00 

angle.  The  extremum  values  of  t^^  differ  by  an  anple  n/p  er-om  t tiose 
of  the  normal  stresses.  For  the  fracture  oalculat.  ions  t lie  stnt.i'  of 
stress  at  0*0  and  it/2  are  Important.  It  in  known  that,  t he  claivape  r.t  res 
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of  the  crystals  is  at  least  twice  the  maxim\im  shear  stress  (Kelly  [?], 
p.  17)*  Of  the  two  state  of  stress  investigated  it  is  found  that 
the  one  at  0=0  plane  is  the  most  important.  In  Fig.  U the  shear  stress 
is  plotted  as  a function  of  p=kr/a  for  a few  hexagonal  crystals  (Zn, 

Cd,  Apatite,  Mg).  The  elastic  constants  of  these  materials  listed  in 
Table  1 are  taken  from  ref.  [8].  The  shear  stress  for  an  isotropic 
crystal  (with  v=0.3)  is  also  plotted  in  Fig.  k.  The  ratio  of  the 
shear  stress  in  hexagonal  crystals  to  that  of  the  isotropic  solids  may 
be  useful  from  the  point  of  view  of  technological  applications.  This 
is  given  in  Fig.  5 for  the  same  crystals.  Finally  we  give  a plot  of 
the  elastic  energy  as  a function  of  the  radius  of  the  cylinder.  Fig.  6. 
Of  course  as  the  radius  of  the  cylinder  R increases  the  elastic  energy 
also  increases.  We  note  however  that  no  singularity  is  present  either 
in  the  stress  field  or  in  the  energy.  The  usual  singxilarities  present 
in  the  classical  elasticity  solutions  however  appear  when  the  lattice 
parameter  a->-0.  This  is  the  classical  continuum  limit.  In  Table  1 we 
give  the  material  moduli,  the  maximum  shear  stress  and  its  radial 
location.  The  ratio  of  the  attenuation  constants  k is  taken,  by  an 
anology  to  the  wave  propagation  [10,  ch.  6],  as 

Table  2 the  energy  ratio  of  the  hexagonal  crystals  to  the  isotropic 
solids  is  listed  for  various  radii  of  the  cylindrical  specimens. 

Finally,  the  maximum  value  of  the  ratio  may  be  used  to  estimate 

the  theoretical  shear  strength  of  the  crystal.  To  compute  this  we  need 
to  estimate  the  attenuation  constant  k in  a(x'-x).  This  fiinction 
decreases  to  its  one  percent  value  at  n lattice  parameter  distance  if 


k = 2.lU6/n 


lit 


As  an  example,  for  zinc  we  have 

(kb/a)"^  t „ = 0.197  10^^  dvn/cm^  . 

r0  max 


If  we  choose  b=a  r.  516)  then 


ai 


t /Cj.  = — — ■ k = 0.050k 

re  max  l.lt  ^^11 


The  maximum  value  that  Kelly  ([7l,  p.  19)  pivec 


to  /C,  , 

re  max  ltii 


0.03lt  . 


Thus  by  comparing  we  obtain 


k = 0.03lt/0.050  = 0.68  , n = 2.lli6/0.68  = 3. lit 


The  maximum  values  of  t „/Ci . are  shown  on  Table  3 for  various 

re  44 

hexagonal  crystals  for  n=2,  n=2.5  and  n=3.0.  These  results  are 
in  the  right  range  as  predicted  by  other  methods  and  atomic 
considerations.  For  example,  in  the  atomic  theory  of  crystal 
lattices  It  Is  known  that  one  must  take  Into  account  the  Interactions 
of  at  least  eight  closest  neighbours  to  obtain  a result  consistent 
with  experiments  [11,  12].  For  hexagonal  crystals  elghtl> neighbours 
corresponds  to  n=2. 

While  clearly  there  Is  one  parameter  (namely  k)  which  need  be 


ii 

T 


estimated,  the  range  of  this  parameter  can  be  ascertained  from  our 
knowledge  in  condensed  matter.  Tlie  detail  shape  of  the  nonlocal 


moduli  seems  to  be  less  effective  so  long  as  it  is  n candidate  to  be 
a distribution.  Tiie  flexibility  in  the  choice  of  a(x'-x)  and  k 
should  be  considered  on  asset  In  the  sense  that  for  noncrystalline 
materials  and  imperfect  crystals  the  attenuation  function  can  be 
adjusted  for  a given  material  once  and  for  all.  Afterward  all 
problems  for  such  a solid  are  reduced  to  boundary-value  problems. 
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Table  1.  Maximum  Shear  Stresses 


Material 


Elastic  Constants  xlO  d 


Mp 

Cd 

Apatite 

Ice(257K) 


2.00  0.617 

6.63  0.O12 


.217 
.3QB  iM 
.237  2.05 

.389  1.55 


6.63  0.O12  0.389  1.55 
0.313  I.0I4  0.396  I.I45 


Table  2.  Strain  Energy  Ratio  1/T. 


p=kr/a 


Material 


Isotropic (v=. 3) 

Isotropic  (ic=l 
v=.3) 


Zn 

Mg 

Cd 

Apatite 

Ice 


0 

5. 

10. 

2.000 

1» . 355 

6.265 

2.001 

U.35I* 

6.263 

0.88q 

2.751 

1+.563 

1.956 

U.20I4 

6.021 

O.8I46 

2.1+56 

3.960 

2.338 

5.272 

7.679 

1.777 

3.807 

5. 1+1+8 

10.121 


iLf  Aft  ^ 9*  ^ n } 

n 

Mtivor  1 fvi 

2.0 

2.5 

8.0 

Zn 

0.05:^ 

0.0l*3 

0.086 

Mg 

0.102 

0.082 

0 . 068 

Cd 

0.071 

0.057 

0.0l»7 

Apat lie 

0.07S 

0.060 

0.050 

Icp 

0.112 

O.O'JO 

0.075 
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FTGllRE  CAPTIONS 


Flf.’rure  I:  (a)  Hexagonal  crystal 

(b)  Straight  edge  dislocation 

Figure  ?:  Cylindrical  coordinates 


t' 


I 


Figure  ?:  Stresses  in  xy  plane  versus  O,  p=P.10. 


Figure  It:  t ./t  versus  p,  9=0.  (v=0.3  for  Isotropic  case) 
ro  o 


Figure  S;  t „/t^,,  , ratio  of  anisotropic  to  Isotropic 
ru  ra 

sliear  stress  at  0=0.  (v=0.3  for  isotropic  case) 


Figure  6:  P/Lb*"  versus  p=kR/a 


Table  1:  Elastic  constants  ami  maxlmiun  shear  stresses 


Table  2:  Strain  energy  ratio 

Table  3:  Shear  stress 
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